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Abstract. We present high resolution (AA ~ 3.7km.s 1 ) 
HST-GHRS observations of the DA white dwarf G191- 
B2B, and derive the interstellar D/H ratio on the line of 
sight. We have observed and analysed simultaneously the 
interstellar lines of Hi, Di, Ni, Oi, Sill and Si III. We de- 
tect three absorbing clouds, and derive a total Hi column 
density N(H l)=2.4±0.1 x 10 18 cm~ 2 , confirming our Cy- 
cle 1 estimate, but in disagreement with other previous 
measurements. 

We derive an average D /H ratio over the three absorb- 
ing clouds N(Di) total /N(Hi) total =1.12±0.08 x 1CT 5 , in 
disagreement with the previously reported value of the 
local D/H as reported by Linsky et al. (1995) toward 
Capella. We re-analyze the GHRS data of the Capella 
line of sight, and confirm their estimate, as we find 
(D/H)ca P eiia = 1.56 ±0.1 x 1CP 5 in the Local Interstellar 
Cloud in which the solar system is embedded. This shows 
that the D/H ratio varies by at least ~ 30% within the 
local interstellar medium. 

Furthermore, the Local Interstellar Cloud is also de- 
tected toward G191-B2B, and we show that the D/H ratio 
in this component, toward G191-B2B, can be made com- 
patible with that derived toward Capella. However, this 
comes at the expense of a much smaller value for the D/H 
ratio as averaged over the other two components, of or- 
der 0.9 x 10 -5 , and in such a way that the D/H ratio as 
averaged over all three components remains at the above 
value, i.e. (D/H) Total = 1-12 x 10~ 5 . 

We thus conclude that, either the D/H ratio varies 
from cloud to cloud, and/or the D/H ratio varies within 
the Local Interstellar Cloud, in which the Sun is embed- 
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ded, although our observations neither prove nor disprove 
this latter possibility. 

Key words: Stars: Individual: G191-B2B - ISM: abun- 
dances - Ultraviolet: ISM 



1. Introduction 

Deuterium is only produced in primordial Big Bang nu- 
cleosynthesis (BBN), and destroyed in stellar interiors 
(Epstein, Lattimer & Schramm 1976). Hence, any abun- 
dance of deuterium measured at any metallicity should 
provide a lower limit to the primordial deuterium abun- 
dance (Reeves et al. 1972). Deuterium is thus a key ele- 
ment in cosmology and in galactic chemical evolution {e.g., 
Steigman, Schramm & Gunn 1977; Audouze & Tinsley 
1976; Vidal-Madjar & Gry 1984; Boesgaard & Steigman 
1985; Olive et al. 1990; Vangioni-Flam & Casse 1995; 
Prantzos 1996, Scully et al. 1997). The primordial abun- 
dance of deuterium is indeed the best probe of the bary- 
onic density parameter of the Universe Qb, and the de- 
crease of its abundance with galactic evolution traces, 
amongst other things, the amount of star formation. 

The first, although indirect, measurement of the deu- 
terium abundance of astrophysical significance was car- 
ried out through 3 He in the solar wind, leading to D/H~ 
2.5 ± 1.0 x 10~ 5 (Geiss & Reeves 1972), a value repre- 
sentative of an epoch 4.5 Gyrs past. The first measure- 
ments of the interstellar D/H ratio, representative of the 
present epoch, were reported shortly thereafter (Rogerson 
& York 1973). Their value of D/H~ 1.4 ± 0.2 x 10" 5 has 
not changed ever since. The most accurate measurement 
of the interstellar D/H ratio was reported by Linsky et 
al. (1993, 1995, hereafter L93, L95) in the direction of 
Capella, using HST-GHRS, D/H~ 1.6 ± 0.1 x 10~ 5 (sta- 
tistical + sytematic). 
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Up to a few years ago, these measurements were used 
to constrain BBN in a direct way. The situation has 
changed, as measurements of the D/H ratio in metal- 
deficient quasars absorbers, at moderate and high redshift, 
have become available (e.g., Carswell et al. 1994; Songaila 
et al. 1994; Tytler, Fan & Buries 1996; Webb et al. 1997; 
Buries & Tytler 1998a,b; see also Buries & Tytler, 1998c 
for a review). However, these observations have not pro- 
vided a single definite value of the primordial D/H ratio. 
At the present time, it is not known whether the higher es- 
timates of the primordial D /H ratio reported are artifacts 
due to the mimicking of the D I line by an H I interloper, 
or whether the ratios reported (low or high) are discre- 
pent due to errors in interpreting the velocity structure, 
or whether substantial fluctuations (a factor ~ 6) actually 
exist. 

On similar grounds, it turns out that determinations 
of the interstellar D/H ratio do not generally agree on a 
single value, even in the very local medium (Vidal-Madjar 
et al. 1986, Murthy et al. 1987, 1990). For instance, D/H< 
10~ 5 is measured toward A Sco (York, 1983), while in the 
opposite direction toward a Aur, D/H=1.65 x 10~ 5 (L93, 
L95). On longer pathlcngths, D/H~ 7. x 10~ 6 is measured 
toward 5 and e Ori (Laurent et al. 1979), and D/H~ 5. x 
1CT 6 toward 9 Car (Allen et al. 1992). Although several 
scenarios have been proposed to explain these putative 
variations {e.g., Vidal-Madjar et al. 1978; Bruston et al. 
1981), the above measurements are still unaccounted for. 
Despite the high quality of the data, it may nevertheless 
be possible that even better data is required to overcome 
systematic effects. 

The only way to derive a reliable estimate of the inter- 
stellar D /H ratio is to observe the atomic transitions of D 
and H in the far-UV, in absorption in the local ISM against 
the background continuum of cool or hot stars. These ob- 
servations have been performed using the Copernicus and 
IUE satellites, and now the Hubble Space Telescope. Both 
types of target stars present pros and cons. 

The main advantage of observing cool stars is that 
they can be selected in the vicinity of the Sun. This re- 
sults in low Hi column densities, and simple lines of sight. 
A difficulty inherent to the cool stars approach is that 
the detailed structure of the line of sight can be obtained 
only through the observation of the Fen and the Mgn 
ions, which are unfortunately not proper tracers of Hi. In 
particular, species like Ni and Oi could not be observed. 
Moreover, the estimate of the H I column density always 
depends strongly on the modeling of the chromospheric 
Lyman a emission line. The Capella target of L93, L95 
is a double system with two cool stars: it provided a line 
of sight with a single absorbing component (see however 
Sec. 5), and allowed a very accurate estimate of the D/H 
ratio, as the emission line could be modeled by observing 
the binary system at different phases. 

Hot stars are unfortunately located further away from 
the Sun, so that one always has to face a high Hi column 



density and often a non-trivial line of sight structure. In 
these cases, Di could not be detected at Lya, and one 
has to observe higher order lines, e.g. Lyy, LyS, Lye. The 
stellar continuum is however smooth at the location of 
the interstellar absorption and, moreover, the Ni triplet 
at 1200 A as well as other Ni lines are available to probe 
the velocity structure of the line of sight. Ni and Oi were 
shown to be reliable tracers of Hi in the ISM (Ferlet 1981; 
York et al. 1983). 

We introduced in Cycle 1 of HST a new type of target, 
white dwarfs, which should solve many of the intrinsic 
difficulties of the problem. Indeed, such targets may be 
chosen in the high temperature range where the depth of 
the photospheric Lyman a line is reduced, so as to obtain a 
significant flux at the bottom of the H I stellar absorption 
line where the interstellar D I and Hi lines appear, as well 
as a smooth stellar continuum. Also, these targets may 
be chosen close to the Sun so that the line of sight is 
not too complex. We observed in HST Cycle 1 the white 
dwarf G191-B2B, at medium resolution AA ~ 18km. s -1 
(Lemoine et al. 1996). 

We proceeded further in this approach with Cycle 5 
high resolution (AA ~ SJkm.s" 1 ) Echelle-A data of Hi, 
Di, as well as Ni and Oi, and Sin and Si in, in order to 
derive an accurate velocity and ionization structure of the 
line of sight. We present these new spectroscopic observa- 
tions of the white dwarf G191-B2B in Section 2. In Section 
3 we present the analysis of these data; the results are fur- 
ther analysed in Section 4. We re-analyze the Capella line 
of sight in Section 5, and summarize our conclusions in 
Section 6. 



2. Observations and Data Reduction. 

2.1. Observations 

Our new observations of the white dwarf G191-B2B were 
performed with GHRS/HST in July 1995 (Cycle 5 Guest 
Observer proposal ID5893). The spectra were acquired in 
the wavelength ranges 1196-1203A, 1200-1207A, 1212- 
1219A and 1299-1307A using the Echelle-A grating. The 
log of these observations is presented in Table [|. 

The Echelle-A grating provides a nominal resolving 
power of ^80 000, or a spectral resolution of 3.7 km.s -1 . 
We used only the Small Science Aperture (SSA), corre- 
sponding to 0.25" on the sky and illuminating 1 diode 
to achieve the best possible resolving power. For further 
details on the instrumentation, see Duncan (1992). 

G191-B2B is a DA spectral type white dwarf located 
40-70 pc away from the Sun, with an effective temperature 
T ofr = 61190-61700K, gravity Log(g)=7.49-7.61 (Finley 
et al. 1997; Vauclair et al. 1997), and magnitude m 1 ,=11.8, 
as determined from optical data fitted with pure hydrogen 
models. It fulfills all the criteria for a good candidate to 
measure the D/H ratio. Apart from Di and Hi at Lyman 
a, we obtained the spectra of the lines Ni(1199A, 1200A, 
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Table 1. List of lines observed in 1995 during Cycle 5 at high spectral resolution AA ~ 3.5km.s 



File 


Element 


Central 


Observing 


Grating 


Exp. 


S/N 






Wave. 


Date 




Time 








(A) 


(1995) 




(s) 




z2q30208m 


Hi, Di 


1215.45 


Jul 26, 22:18 


ECH-A 


6528 


21 


z2q3020ct 


Hi, Di 


1215.45 


Jul 27, 00:59 


ECH-A 


6528 


21 


z2q3020gt 


Hi, Di 


1215.45 


Jul 27, 03:39 


ECH-A 


6528 


21 


z2q30108t 


Ni 


1199.94 


Jul 26, 13:30 


ECH-A 


4896 


25 


z2q3010dt 


Ni, Sim 


1203.57 


Jul 26, 18:27 


ECH-A 


3917 


19 


z2q3010cm 


Oi, Sin 


1303.01 


Jul 26, 16:31 


ECH-A 


4570 


26 



1201A), Sim(1206A), Ol(1302A) and Sill(1304A) using 
Echelle-A. The observation at high spectral resolution of 
O I and N I should allow us to resolve possible different H I 
absorbing clouds on the line of sight, while the observa- 
tion of Si II and Si in should allow the detection of possible 
H II gas. From the comparison of the b values of different 
atomic weights, an estimate of the temperature and turbu- 
lent velocity in the different clouds is also possible, while 
the analysis of the Sin and Sim lines should determine 
the ionization structure of the line of sight. In particular, 
the observed lines of Ni, Sin and Si III are not saturated, 
and should therefore provide accurate estimates of the col- 
umn densities. The Oi 1302A line is slightly saturated in 
this low Hi column density environment, and should thus 
offer a reasonable estimate of the column densities, as well 
as of the b- values. At the end, we expect to obtain a suffi- 
cient number of constraints on all parameters involved in 
the analysis of the Lyman a lines of D I and H I, but the 
D/H ratios of the different components. 

2.2. Data reduction. 

Our data were reduced with the Image Reduction and 
Analysis Facility (IRAF) software, using the STSDAS 
package. During the observations, we used the FP-SPLIT 
mode which splits the total exposure time into successive 
cycles of 4 sub-exposures, each corresponding to a slightly 
different projection of the spectrum on the photocathode. 
We used the "quarter stepping" mode, which provides a 
sample of 4 pixels per resolution element. This oversam- 
ples the spectrum (since for instance the SSA does not ful- 
fill the Nyquist sampling criterion), and allows to correct 
for the granularity of the photocathode. Indeed, the effect 
of the photocathode on each diode being the same for the 
four sub-exposures, it is possible to evaluate this granu- 
larity from the comparison of the 4 sub-exposures where 
a constant granularity effect mixes with a non-constant 
photon statistical noise. 

We found that the standard method for correcting for 
the granularity, which is available in the IRAF-STSDAS 
package, was not efficient as a result of the impact of the 
geomagnetic field on the photocathode. The geomagnetic 
field deflects the incident electron beam on the diodes in 
a way depending on the location of the HST on its orbit, 



thus making the granularity pattern non-constant from 
sub-exposure to sub-exposure. As a result, the standard 
iterative method leads to the appearance of noisy features 
whose intensity increases with the number of iterations. 
We therefore developed a different procedure based on the 
use of simple statistical filters. 




1212 1214 1216 1218 1220 
WAVELENGTH 

Fig. 1. Data taken with the Echelle-A in the Lyman a 
region. Some perturbations due to the photocathode in- 
homogeneities can be seen near 1213.4 and 1217.5 A. Note 
the deuterium absorption at 1215. 3A. The central peak at 
the bottom of the saturated H I Lyman a line corresponds 
to the earth geocoronal emission. Note also that the bot- 
tom of the saturated Lyman a line is not at exact zero 
level. 

We first cross-correlate the different spectra corre- 
sponding to the same spectral ranges, i.e. the different 
sub-exposures of a same exposure, and shift them accord- 
ingly so that real absorption or emission features corre- 
spond to the same pixel values in each spectrum. We next 
compute the median and the standard deviation with re- 
spect to that median at each pixel, our statistical sample 
being the set of intensities of the different spectra at that 
pixel. Finally, we compute the average intensity value at 
each pixel, rejecting the intensities of the different spectra 
deviating from the median by more than, say, 3 standard 
deviations. We found that this simple treatment resulted 
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in good signal-to-noise ratios (Fig. ^) while avoiding any 
residual ghost feature (see Bertin et al. 1995 for more de- 
tails). 

We chose to rebin our spectra to two pixels per diode in 
all cases but O I, i.e. two pixels per resolution element ap- 
proximately, which is the standard sampling mode, using 
the IRAF-STSDAS procedure. In the case of O I the very 
sharp edges of the line contain most of the information, 
and we thus decided to keep four pixels per resolution ele- 
ment so as to take advantage of the oversampling available 
in these observations. 

In the case of Lyman a, the spectral signature is not as 
sharp, and, moreover, the S/N ratio is significantly lower 
at the bottom of the line, and we thus developed a sec- 
ond approach to check the consistency of our procedure. 
We re-align the sub-exposures using the deuterium line, 
which is clearly detected in each of the sub-exposure, as a 
reference. This procedure means we are unable to correct 
for the photocathode granularity. Nonetheless, it appears 
that no such defect is present at Lyman a. In effect, such 
defects are of order ~10%, and are thus easily detected 
when adding the sub-exposures corresponding to a priori 
the same spectral instrument shift, thereby improving the 
S /N ratio and building four different shifted spectra where 
the photocathode defects appear at fixed positions. This 
process allowed us to show that the amplitude of photo- 
cathode defects in the Lyman a central region does not ex- 
ceed a few percent. One clearly sees such defects in Fig. 
in the far wings of the Lyman a profile, around 1213.4 
and 1217. 5A (these defects were not perfectly corrected 
through the first approach). Thus the final Lyman a line 
profile is well represented in the region from 1213.6A to 
1217. 3A, and the region that we use in our data analy- 
sis is even further reduced to 1214. 5-1217. OA, where the 
spectral information lies. 

As a second consistency check, we measure the width 
of the geocoronal emission line. Spurious shifts in the co- 
adding procedure would degrade the spectral resolution. 
The geocoronal emission line should be a very narrow fea- 
ture, only slightly broader than the nominal instrument 
spectral resolution. We find, fitting the emission line with 
a gaussian, an equivalent Doppler width b ~ 5.7km.s~ , 
corresponding to a temperature T ~ 1950K, after decon- 
volution from the instrumental profile. At this point, one 
should note that the geocorona fills up the SSA, so that its 
width results from the convolution of the intrinsic width 
with the spectrograph instrumental profile, and with the 
profile of a point source broadened by the projection of 
the SSA on the spectral scale; this latter effect amounts 
to broadening the SSA instrumental profile by 4 pixels. 
The above measured width thus corresponds to an in- 
trinsic temperature T ~ 600K, which is a very plausible 
average value for the Earth exosphere temperature. We 
note that this consistency test is all the more significant 
as the geocoronal emission line only shows up in the fi- 
nal co-added sub-exposures. Therefore, the procedure used 



for co-addition produces a sharp line profile in which the 
Echelle-A spectral resolution is kept at its highest nominal 
value. 
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Fig. 2. The normalized Nl triplet data used in the fitting 
procedure. 
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Fig. 3. Same as Fig. [| The Oi line. 



Each exposure was preceded by a special platinum 
lamp calibration which allowed us to calculate the resid- 
ual shift left by the standard CAL-HRS calibration proce- 
dure. Together with the oversampling mode, this correc- 
tion allowed us to reach an absolute calibration accuracy 
of il.5km.s~ 1 in radial velocity. Hereafter, the radial ve- 
locities will be given in the heliocentric rest frame. 

The different spectral lines were normalized to a con- 
tinuum of unity by a polynomial whose degree, typically 
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Fig. 4. Same as Fig. | The Sin line. 
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Fig. 5. Same as Fig. g. The Si III line. Note that the broad 
line on the red side is of photospheric origin (see text). 



2 — 7, was chosen according to a cross-validation statisti- 
cal procedure. The normalized spectra are shown in Fig. ^], 
Fig. [| Fig. [| and Fig. ||. Since the continuum level is de- 
termined to high precision, we do not need to take into 
account any continuum fitting parameters in the line fit- 
ting. The Lyman a line was not normalized in this way, as 
the stellar continuum is a Lorentzian profile, with a deep 
core; the stellar continuum at Lyman a is described at 
length later. 

Finally, in the case of the Si ill line (1206A), we de- 
cided to keep in the spectrum the nearby signature of the 
photospheric Si in stellar line on the red wing of the inter- 
stellar Sim line (see Fig. ||). This photospheric line was 
thus fitted along with the interstellar features as an ad- 
ditional component. This gave access to an independent 




10 20 
Heliocentric Velocity [km/s] 

Fig. 6. The absorption features seen in different species, 
as indicated on the left hand side of the figure, superposed 
in velocity space. Two obvious components are seen in Si n 
and Ni. A third component in between the other two is 
detected in the Si in and Ol lines. The two components 
at the shortest wavelengths, seen in Si ill, are Hn regions, 
while the component at ~20 km s" 1 is the Local Cloud, an 
Hi region. It is obvious that deducing velocity structure 
of lines of sight for D/H evaluations from ionized species 
could be extremely misleading. 



evaluation of the intrinsic shift of the white dwarf photo- 
spheric features, a parameter that is also involved in our 
final analysis of the deuterium and hydrogen Lyman a 
profiles. 

2.3. The number of absorption components along the line 
of sight 

We first evaluate the number of components present on the 
line of sight. This number represents a minimum number 
of components since blends are always possible. 

Figure || illustrates all species (excluding H I and D i) 
observed with the Echelle-A spectrograph. It is clear, from 
the N I and Si II profiles, that at least two clearly resolved 
components are present. In the Si III line, the strongest 
component falls in between the other two seen in N I and 
Si II. This cannot be due to an incorrect absolute cali- 
bration of one of the spectral ranges, as we observed the 
Sim line at 1206A together with the Ni line at 120lA in 
the same exposure. The relative shift between these two 
species, from one end of the detector to the other, should 
certainly be less than 0.5 km.s^ 1 . This Sim component is 
thus real, and obviously more ionized than the other two. 

A broad component is also detected in Si in at the edge 
of the observed spectral range. When fitted with a single 
absorption feature, assuming that the broadening is ther- 
mal, the corresponding temperature is of order 30000 K, 
i.e. much broader than clouds usually seen in the local 
ISM. Because its heliocentric velocity is ~ 29km. s _1 , it is 
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reasonable to assume that this line is of photospheric ori- 
gin. Indeed, as shown in Vidal-Madjar et al. (1994), pho- 
tospheric features in G191-B2B are relatively broad and 
their observed heliocentric velocities are ~ 35 ± 4km. s -1 . 
Predictions of the strength of the Si in photospheric line 
are also in agreement with the observations. We thus do 
not include this feature in our study. 

The O I profile confirms, through profile fitting, the 
need for a third insterstellar component, located in be- 
tween the two components seen in N I and Si II. The O I sig- 
nature of this third component is mainly the flatness of the 
absorption in the saturated core. We also note that there 
is no component around 29 kru.s - in the other lines, con- 
sistent with the photospheric interpretation of the broad 
Si in feature. 

2-4- The residual flux at the bottom of the spectral lines. 

In order to properly evaluate the strength of the different 
ISM components contributing to the observed absorption, 
it is extremely important to evaluate the level of the zero 
flux which may be slightly erroneous in an echelle spectro- 
graph, due to the diffuse light produced from the adjacent 
spectral orders. The case of Lyman a is obvious as the 
broad saturated core allows us to accurately determine 
for the zero level. We corrected the Lyman a region for 
this zero level by subtracting the flux in the Lyman a core 
from the whole line profile. Correcting for the zero level in 
the GHRS data reduction pipeline is indeed difficult for 
other species, as it depends on instrumental parameters 
as well as on the (unknown) spectrum shape in adjacent 
orders. Thus, for species other than Hi, we determined the 
zero level as the best fit value, i.e. the one corresponding 
to a minimum x 2 , of all Oi, Ni, Sin, and Sim lines fitted 
simultaneously. In each x 2 calculation as a function of the 
zero level, the same zero level value was adopted for the 
O I and Si n lines, and for the N I and Si in lines, since 
these pairs were observed in the same sub-exposures. In 
each case, only the Oi and Ni lines contributed signifi- 
cantly to determining the zero level values, as the lines 
of Si II and Si in are too weak. The N I lines are sensitive 
to the zero level as the strongest line is almost saturated, 
and because these lines form a triplet and their measured 
relative strengths depend on the zero level. For our data 
normalized to unity, the value and the error bar of the N I 
and Si in zero levels is: 0.038 ± 0.02, i.e. consistent with 
zero level within the accuracy of our measurement. Since 
the O I line is saturated, the zero level position is sharply 
defined. The measured value of the zero level, for both 
Oi and Sin lines, is: 0.032 ± 0.004. Theses values will be 
used from now on in our analysis and their impact on the 
evaluated parameters will be considered as minor in the 
forthcoming study and checked in few specific cases. 

3. The total Hi and Di content 



3.1. The stellar Lyman Alpha profile 

It is well known that the limiting factor of an accurate 
D/H ratio measurement is generally the measurement of 
the Hi column density itself, as the Di line at Lyman a 
is not saturated for N(Hi)^ 10 18 cm -2 . In the present 
case, N(Hi) can be measured with good accuracy from 
the damping wings of the H I Lyman a profile. How- 
ever, we first need to model the stellar continuum, as the 
Lorentzian shape of the intrinsic Lyman a photospheric 
line may mimic interstellar damping wings, and affect the 
N(Hi) measurement. We thus calculated Lyman a pro- 
files from LTE atmosphere models using the appropriate 
stellar parameters. Another LTE profile, as well as NLTE 
calculations, both using codes completely different from 
ours, were kindly provided by Stefan Drcizler. A qualita- 
tive difference between these profiles is that LTE profiles 
are not as sharp and deep as NLTE profiles. 

In our analysis, the different possible continuum were 
considered (see Fig. in order to evaluate the possible 
systematic error one could make on the total Hi evalua- 
tion. 



G191— B2B Lyman Alpha 
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Fig. 7. The different calculated Lyman a photospheric 
continua assumed in our analysis: from top to bottom, a 
simple 2 nd order polynomial, a first LTE calculation (DK), 
a second LTE calculation (S. Dreizler, private communica- 
tion) and an NLTE calculation (S. Dreizler, private com- 
munication). These three latter profiles were shifted by 
29 km.s -1 to match the measured radial velocity of G191- 
B2B. 

The heliocentric radial velocity of G191-B2B photo- 
spheric lines (including gravitational redshift) is known to 
be 24 ± 5 km.s -1 from ground based observations (Reid 
and Wegner 1988). From our Cycle 1 observations, we 
measured a differential velocity ~ 35±4km.s~ 1 compared 
to the nearby interstellar lines (Vidal-Madjar et al. 1994); 
this shift was however evaluated from medium resolution 
data in which the ISM components were not resolved. 



Vidal-Madjar et ai: Detection of local ISM D/H variations. 



7 



In the present high resolution Echelle-A data, the photo- 
spheric Si in line falls at a radial velocity 29 ± l.Skm.s - . 
This line is however located at the end of our spectral 
range, and the profile shape may be partially distorted, 
inducing a shift in the measured velocity. We also tenta- 
tively identified six Fe V and six Ni v photospheric lines 
in the 1300-1306A region, whose average location gives a 
radial velocity 24.8 ± 4.0km.s" 1 . All these measurements 
are in good agreement, and their weighted least square av- 
erage results in a radial velocity (including gravitational 
redshift): 28.9 ± l.Skm.s" 1 . With a gravitational redshift 
of order ~ 25km. s" 1 (Finley et al. 1997), the radial veloc- 
ity of G191-B2B itself is thus of order ~ 5km. s" 1 atypical 
value for an object in the local galactic environment. 

With this information, we produced different data sets 
for the Lyman a line, each corresponding to the normal- 
ization of the observed data by a theoretical stellar contin- 
uum. From now on, we fix the radial velocity of the stellar 
profile at 29km. s" 1 . Finally, we defined an additional "Ly- 
man a photospheric profile" with a second order polyno- 
mial, in order to circumvent any assumption made about 
the shape of the profile itself. This polynomial continuum 
is derived by fitting to the far wings of the interstellar line 
in the ranges 1213.6-1213.9 A and 1217.0-1217.3 A (see 
Fig. Q). In fact, we found that the best solution (minimum 
X 2 ) was obtained for this second order polynomial further 
corrected by a fourth order polynomial in the core of the 
line, the latter being fitted together with the interstellar 
Lyman a line. 

We now have four different profiles, hence four normal- 
ized data sets. From these, we discard the stellar profile 
corresponding to the LTE calculation of S. Dreizler, since, 
as shown in Fig. ^, this profile is intermediate between our 
LTE calculation, and the NLTE calculation of S. Dreizler. 
We thus keep these two latter extreme case, plus the sec- 
ond order polynomial. Hereafter, we respectively denote 
these continua as LTE29, NLTE29, and 2 nd Polynomial. 

3.2. The total Hi and Di column densities 

We now evaluate the total column densities of Hi and 
D I along the line of sight. We make different assumptions 
about the stellar continuum as well as the number of com- 
ponents along the line of sight, in order to evaluate their 
impact on the estimate of the column densities. Namely, 
for each of the above three stellar continua, we assume 
a further correction to the continuum, given either by a 
first order polynomial, or a fourth order polynomial; in 
each case, these further corrections are fitted simultane- 
ously with the physical parameters of the line of sight. 
As to the number of components on the line of sight, we 
performed the calculation for one, two, three, and four in- 
terstellar components. We certainly know from the other 
elements that there are at least three absorbing compo- 
nents; what we want to show here, however, is that the 
evaluation of the total column densities, and hence the 



Table 2. Reduced x 2 (with 349 degrees of freedom) and 
several parameters corresponding to the Lyman a fits with 
the 3 different assumed stellar profiles: 2 nd order polyno- 
mial, LTE29 and NLTE29 (see text). Only the H I and D I 
lines were fitted, with all parameters free (see text), and 
only one component on the line of sight. The order of the 
polynomial adjusted on top of the reconstructed stellar 
continuum is indicated, either l st order or 4 th order. 



1 st order polynomial fit 


Stellar cont. 2 nc, Pol. 


LTE29 


NLTE29 


X" 1-91 


1.51 


2.52 


N(Hi) To t (cm" 2 ) 3.00 x 10 18 


2.21 x 10 18 


1.80 x 10 18 


N(Dl) To t (cm" 2 ) 2.43 x 10 13 


2.51 x 10 13 


2.66 x 10 13 


(D/H) Tot (xlO 5 ) 0.81 


1.14 


1.48 


4 th order polynomial fit 


Stellar cont. 2 nd Pol. 


LTE29 


NLTE29 


X J 1-13 


1.51 


2.02 


N(Hl) To t (cm" 2 ) 2.32 x 10 1S 
N(Di) T ot (cm" 2 ) 2.60 x 10 13 


2.21 x 10 18 


2.00 x 10 18 


2.66 x 10 13 


2.66 x 10 13 


(D/H) T ot (xlO 5 ) 1.12 


1.20 


1.33 



Table 3. Same as Table ^ but with three components 
on the line of sight. All three components are required to 
present a unique D/H value. Here, the total D/H ratio 
represents the ratio of the total D I column density to the 
total H I column density. 



1 st order polynomial fit 


Stellar cont. 


2 nd Pol. 


LTE29 


NLTE29 


x 2 


1.83 


1.33 


2.05 


N(Hl) T ot (cm" 2 ) 


2.92 x 10 18 


2.12 x 10 18 


1.70 x 10 18 


N(Di) To t (cm" 2 ) 


2.58 x 10 13 


2.70 x 10 13 


2.93 x 10 13 


(D/H) Tot (xlO 5 ) 


0.88 


1.27 


1.72 


4 th order polynomial fit 


Stellar cont. 


2 nd Pol. 


LTE29 


NLTE29 


x 2 


1.03 


1.32 


1.71 


N(Hl) To t (cm" 2 ) 


2.38 x 10 18 


2.01 x 10 18 


1.92 x 10 18 


N(Di)tm (cm" 2 ) 


2.73 x 10 13 


2.73 x 10 13 


2.79 x 10 13 


(D/H) Tot (xlO 5 ) 


1.15 


1.30 


1.45 



average D/H ratio, are not sensitive to the number of com- 
ponents. Finally, in order to further show the strength of 
this prediction, we do not take into account any physical 
constraint coming from other species, i.e. Oi, Ni, Sin, 
or Si in. When those elements are included in a global fit 
of the line of sight, this will narrow down the margin of 
freedom for the D/H ratios. 

The results of some of the parameters of these fits 
are presented in Table |2| for a one cloud solution, and 
in Table |^ for three absorbing components. The results 
are similar for a two component structure of the line of 
sight. Several comments can be made at this stage. We 
note that the theoretical continua have very well defined 
shapes, with abrupt features in their profile at various 
points near the core. When the interstellar absorption is 
superimposed, mismatches could be introduced in the re- 
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Fig. 8. The different fits are presented for the Lyman a line when several stellar continua are assumed. The upper 
panels correspond to a I s * order polynomial fit on top of the assumed stellar continuum, while the lower panels 
correspond to a 4 th order polynomial fit. From left to right the pannels correspond to the three stellar continua 
presented here, i.e. 2 nd Polynomial, LTE29 and NLTE29 (see text). The central panels show that the LTE29 profile 
is probably a good approximation to the real stellar photospheric profile, since the 4 th order polynomial fit does not 
improve the fit with respect to the I s * order polynomial. On the contrary, the 2 nd Polynomial is not deep enough in 
the core, while the NLTE29 is too deep. Once corrected for their erroneous shape, the solutions found in terms of total 
Hi and total D/H are similar (see Table ||). The x 2 values favor the 2 nd x 4 th stellar continuum (see text). 



suiting final profile. However, as shown in Fig. ^[ the 4 th 
order correction deepens the 2 nd order polynomial in the 
core of the line to mimic a Lorentzian profile, but hardly 
modifies the LTE29 continuum. This means that the A\ 2 , 
of order 100, in favor of the 2 nd order continuum, is mostly 
intrinsic to the LTE29 profile. 

The above fits also reveal that LTE continua are gen- 
erally favored (in terms of \ 2 ) when compared to NLTE 
continua. This is illustrated by the fact that a 4 th order 
correction plus the NLTE profile tends to the LTE contin- 
uum (see Fig. ^). One could be surprised by this result, 
since NLTE model atmospheres should be more appro- 
priate than LTE calculations at the high temperature of 
G191-B2B. That the LTE profiles seem to fit better may 
partly be due to the fact that we compare pure hydrogen 



models only. G191-B2B is well-known to contain signifi- 
cant amounts of heavy metals, which are visible in high- 
resolution UV spectra and which even completely domi- 
nate the EUV spectrum (e.g., Wolff et al. 1998). Taking 
into account the blanketing effect of these metals leads to 
a decrease of the effective temperature by about 5000 K, 
and may also lead to a better agreement of NLTE profiles 
with the observations. These details are, however, incon- 
sequential for the present study. 

In summary, these x 2 analyses give two main results: 
the 2 nd x 4 th order polynomial is largely favored over the 
theoretical continua; it is also supported by the absolute 
X 2 value, with a reduced x 2 close to unity. We can thus 
conclude that this continuum offers a very good approxi- 
mation to the actual stellar continuum of G191-B2B, and 
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we adopt this for our determination of the (D/H) ratio. 
When we include all other species in the profile fitting, 
we obtain the same result, namely that the 2 nd x 4 th or- 
der polynomial continuum is very largely favored over the 
other theoretical continua. Although all subsequent pro- 
file fittings were also carried out for these continua, we 
will not refer to them any more. Finally, we also checked 
that the presence of more than 3 components on the line 
of sight does not change these results. Indeed, repeating 
the same procedure with 4 components led to results that 
are nearly identical to those of the 3 cloud solution, in 
terms of total Hi and Di column densities. 

Given these considerations, we study the relative vari- 
ation of the x 2 as a function of the total D/H ratio, in 
order to derive confidence intervals on this quantity, ac- 
cording to the so-called A% 2 method. Strictly speaking, 
this method is exact for gaussian quantities, and, more 
generally, correct up to second order in the variation of 
the parameters. The 2a confidence level around a given 
parameter should then be reached when the % 2 of the fit 
with the parameter fixed to this ±2cr value becomes equal 
to \q + 4, where Xo is the total best-fit x 2 with the pa- 
rameter left as free. The 3a level is obtained for Xo + 9> 
etc... However, this assumes that the quantities are gaus- 
sian, that the data points (pixel values) are uncorrelated, 
and that their error bars are accurately known. We prefer 
to relax these assumptions, and set the threshold for an 
"effective" 2a level at Xo + 10- This provides a very conser- 
vative and very safe range of errors. In particular, we note 
that this procedure has the great advantage of including 
in the error bars so-derived all errors on other parameters 
that would have to be propagated and projected on the pa- 
rameter under study. In effect, when looking for a x 2 value 
for a given value of the parameter, away from the best fit 
value, the profile fitting tends to accomodate the values 
of the other parameters so as to obtain the smallest x 2 
hence the smallest A% 2 ■ This is especially important when 
different parameters are strongly correlated, e.g. column 
densities anti-correlated with broadening parameters, or 
column densities correlated with other column densities. 
However, we note that, strictly speaking, this method can 
only be applied to statistical errors, and not to systematic 
errors, which tend to add linearly, rather than quadrati- 
cally. Since this method is CPU time expensive, we only 
use to determine error bars on D/H ratios. For other quan- 
tities, that do not constitute our primary concern, we will 
quote error bars estimates obtained through trial and er- 
ror. 

Here, from the variations of the % 2 as a function of the 
D/H ratio, we obtain, for three components, and a unique 
D/H ratio: 

(D/H) To tai = 1.15 ±0.1 x 1(T 5 . 
In the above, uncertainties arising from the choice of 
the stellar continuum are included in the quoted error 
bars, because of the simultaneous fitted 4 th order correc- 
tion to the continuum. However, systematic errors, that 



could arise from uncertain calibration, or physics left out 
of our model, are not included. 

4. The component by component analysis 

4-.1. The components velocity separation and column 
density ratios 

No information on the velocity structure of the line of sight 
can be obtained from Lyman a as the components have 
too large intrinsic widths in D I and H I to be resolved. 
Therefore, we now turn to a detailed analysis of the other 
species, Oi, Ni, Sill, and Si in, to extract the required 
information. 

In order to do so, we perform a fit of all lines of Oi, 
Ni, Sin, and Si III, simultaneously, assuming three com- 
ponents on the line of sight. This means, more precisely, 
that we define one velocity, one temperature, and one tur- 
bulent broadening for each cloud; the column densities of 
the various species are also free parameters. The absorp- 
tion profile corresponding to these physical parameters is 
then computed for each absorber and for each absorption 
line, and the total profile is compared to the observations. 
We then minimize the sum of the total x 2 of each absorp- 
tion line. The best fit values of the physical parameters 
is presented in Table |[ In this table, most of the error 
bars are tentative, obtained by trial and error, i.e. they 
represent a rough estimate of the actual error bar; this is 
however sufficient for our present purpose. 

Among these physical parameters, the radial veloci- 
ties and the column densities are the best determined, as 
would be expected. The temperature and turbulent broad- 
ening are indeed slightly degenerate as the mass separa- 
tion between these species is insufficient. We note that the 
radial velocities are determined to a high degree of con- 
fidence in the spectral regions where they are resolved, 
e.g. Si III for component 2, and Ni, O I, and Sin for com- 
ponents 1 and 3, due to the large number of sampling 
points on the profile. (We recall that the components are 
numbered from 1 to 3 in order of increasing radial ve- 
locity). The error on the estimate of their radial veloc- 
ity is dominated, in a given spectral region, by the cal- 
ibration accuracy ±1.5 km.s -1 . When averaged over all 
spectral regions, the accuracy should thus be of order 
±0.8km.s _1 . The relative shifts in radial velocity between 
the various components are not affected by the absolute 
calibration accuracy; their error bars are however of or- 
der ±0.5km.s~ 1 but will be greatly improved later. We ob- 
tained these error bars by performing fits to all species 
but one, each species being excluded in turn, and averag- 
ing the resulting radial velocities and their shifts. 

It is especially important to note the radial velocity of 
the third component, V = 20.3 ± 0.8km. s _1 . This compo- 
nent is identified with the Local Interstellar Cloud (LIC) 
in which the Sun is embedded (e.g., Lallement & Bcrtin 
1992; Bertin et al. 1995). This cloud has been detected in 
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Table 4. The ISM components characteristics in the direction of G191-B2B as evaluated from the Ol, Nl, Sill and 
Si in lines (from the blue to the red the components are noted 1, 2 and 3; component 3 is identified with the LIC and 
noted as such). The error bars are estimated from the different fits produced, i.e. fits with all the lines and with all 
the lines but one included. Note that the error bars on components 1 and 2 are in some cases larger simply because 
these two components are blended. As well, note that the values given here will change slightly when the H I and D I 
lines are included in the profile fitting, see subsequent tables. 





Component 1 


Component 2 


LIC 


V (km.s" 1 ) 


8.2 ±0.8 


13.2 ± 0.8 


20.3 ±0.8 


T(K) 


8000. ± 4000. 


8000. ± 6000. 


7000. ± 4000. 


a (km.s 4 ) 


1.5 ±1.5 


2.5 ± 1.5 


1.2 ± 1.2 


N(Oi) (cm- 2 ) 


2.8 ± 1.0 x 10 14 


1.5 ± 1.0 x 10 14 


3.1 ±0.5 x 10 14 


N(Nl) (cm- 2 ) 


1.2 ±0.3 x 10 13 


3.5 ± 1.0 x 10 12 


6.7 ±0.4 x 10 13 


N(Siii) (cm~ 2 ) 


1.9 ±0.3 x 10 13 


8.4 ±4.0 x 10 12 


1.1 ±0.2 x 10 13 


N(Sini) (cm" 2 ) 


4.0 ±3.0 x 10 11 


2.2 ±0.5 x 10 12 


< 10 11 



the direction of nearly all nearby stars, and its velocity 
vector relative to the Solar System was found to be: 

Vlic = 25.7km.s~ 1 

1//(LIC)= 186.1° 

b//(LIC)= -16.4°. 

This component is actually the only absorbing cloud 
detected toward Capella, whose line of sight is only sepa- 
rated by ~ 7° from the G191-B2B line of sight. The pro- 
jections of the LIC velocity on the Capella and G191- 
B2B lines of sight are respectively 21.96km.s~ , and 
20. 26km. s -1 . We therefore obtain a radial velocity in ex- 
cellent agreement with the prediction. This value is of im- 
portance as it allows us to check the consistency of our 
approach: we expect to detect this component, and we 
know its predicted velocity within the calibration internal 
error ±0.8km.s~ 1 . Taking into account the appearance of 
the LIC in four spectral domains (Di, Ol, Nl, and Sin), 
we conclude that we actually well detected and correctly 
identified it with component 3. From now on, we thus refer 
to this component as the LIC. 

Finally, we note that the column density ratios be- 
tween the components are relatively close to 1. This is 
illustrated in the following Table (from the blue to the 
red the components are noted 1, 2 and LIC; component 1 
is used as the reference): 



Components N ratios 


2/1 


LIC/1 


Ol 


0.5 


1.1 


Nl 


0.3 


5.6 


Sin 


0.4 


0.6 


Sim 


5.5 


< 0.02 



This also provides another consistency check. Namely, 
when a solution is obtained for H I and D I, the column 
density ratios between the different components should 
remain close to the values obtained for Ol, and/or Nl. 
As far as Nl, Sin, and Si III are concerned, the error on 
these ratios, and in particular the LIC to 1 ratio, is ex- 
pected to be relatively small, 55 40%, as the lines are not 



saturated. Notably, this implies that the various compo- 
nents have very different ionization structures. Obviously, 
components 1 and 2 are more ionized than the LIC. In- 
terestingly enough, O I and N I do not appear to behave 
similarly; this trend will be confirmed below, where a more 
detailed analysis is performed including H I and D I. 

4-2. The Hi and Di content of each component 

We now turn to the final analysis, i.e. all observed lines 
fitted simultaneously. We recall that our fitting procedure 
works in terms of absorbing components, so that a unique 
set of radial velocity, temperature and turbulent broaden- 
ing is used to calculate all absorption profiles of all lines 
for a given component. These physical parameters are thus 
fitted simultaneously over the lines of D I, Hi, O I, N I, Si n, 
and Si in; we assume the presence of three components on 
the line of sight. The total number of degrees of freedom 
available is ~ 800. 

The best fit, without including any extra constraint 
on the parameter space, is shown for the Lyman a line in 
Fig. [|, and is enlarged around the D I line in Fig. [H]. This 
best fit is obtained for a 2 nd x 4 th order continuum, and 
the reduced \ 2 is X 2 — 1-12 for 795 degrees of freedom. 
We note that the physical parameters obtained, given in 
Table |, are in good agreement with those obtained from 
the fits of the Ol, Nl, Sin, and Si in lines only. Moreover, 
the individual x 2 corresponding to each spectral region, 
derived here from a simultaneous fit of all lines, are close 
to those obtained for individual fits of each spectral re- 
gion. This supports the idea that the structure of the line 
of sight is well determined. Finally, we find column den- 
sity ratios between components, in H I and in D I, that lie 
relatively close to those obtained for O I and N I in the pre- 
vious fit of O I, N I, Si II, and Si III, as well as in this global 
fit. That our results indicate that Hi, Di, Ol, and Nl 
behave in a similar way, represents a strong consistency 
check. Indeed, it is well known that these four elements 
have very similar ionization properties, and, in particular, 
that they are locked by charge exchange in a typical ISM 
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Table 5. Reduced % 2 , A% 2 (total) and several parameters corresponding to the fits made with the different constraints 
assumed. All observed lines are included in the profile fitting, and all parameters are free, except otherwise noted in 
boldface. The first column gives the best-fit all parameters free result; the subsequent columns give results for various 
constraints: respectively, for a unique D/H ratio between the three components, for a unique N(Oi)/N(Hi) ratio, for 
a unique N(Ni)/N(Hi) ratio, for a D/H ratio in the LIC corresponding to that found by L93, L95, and, finally, for the 
LIC D/H ratio, temperature and turbulent broadening corresponding to those found by L93, L95, toward Capella. In 
all cases, the stellar continuum is a 2 nd x 4 th order polynomial. 
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N(Di) L ic(x10- 13 ) (cm" 2 ) 


1.10 


0.940 


1.10 


1.42 


1.10 
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(D/H) L ic (xlO 5 ) 


1.10 


1.09 


0.94 


0.86 


1.60 


1.60 


AV2-1-AV3-1 (km.s" 1 ) 


5.03-12.18 


5.07-12.14 


5.16-12.24 


5.03-12.05 


5.14-12.04 


5.11-12.05 


Hi ratio 2/1-3/1 


3.48-3.22 


4.58-3.14 


0.59-1.61 


0.32-5.11 


6.91-3.14 


6.76-3.52 


Dl ratio 2/1-3/1 


1.96-2.06 


4.58-3.14 


2.31-2.31 


1.39-2.16 


2.05-2.10 


2.20-2.47 


Ol ratio 2/1-3/1 


0.57-1.12 


0.56-1.10 


0.59-1.61 


0.56-1.26 


0.49-1.10 


0.52-0.91 


Nl ratio 2/1-3/1 


0.41-5.61 


0.40-5.61 
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0.32-5.11 


0.34-5.37 
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0.44-0.52 
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4.99 


5.35 


4.99 


4.25 


4.77 



environment. Provided that these elements have a unique 
abundance on the line of sight, one should therefore expect 
the abundance ratios of their neutrals to be unique. 

We now go one step further with this idea, and test 
whether these elements can actually afford a unique neu- 
tral abundance on the line of sight. We thus require, 
in turn, the ratios N(D i)/N (H 1), N(Oi)/N(Hi), and 
N(Ni)/N(Hi) to be unique (i.e., equal from component to 
component), and the unique abundance kept as a free pa- 
rameter, all other parameters as equally free as above, and 
we compare the % 2 of these three solutions with the above 
best-fit x 2 - The results are given in Table |. Obviously, 
Dl can be considered as a perfect tracer of Hi on this 
line of sight, as the relative x 2 difference with the above 
best fit is reasonably small, A\ 2 = 3.9. Two remarks are 
in order here. This result means that, within the quality 
of our data set, we do not detect variations of the D/H 



ratio, from cloud to cloud; it does not mean either that 
we rule out such variations. Furthermore, that Dl seems 
to trace Hi holds here because we only consider one line 
of sight. When we compare the value for the D/H ratio 
measured here toward G191-B2B, with that measured to- 
ward Capella, in the next section, we will find that they 
disagree, and thus argue that the D/H ratio varies within 
the local ISM. 

For 1, the situation is not as clear, as the x 2 difference 
is larger, A% 2 = 14.6. However, we note that in this profile 
fitting, the instrumental zero-level of the O I line was kept 
fixed (for numerical reasons) . The 1 line is saturated, 
and the column density ratios between components in 1 
thus depend on the value of the zero-level. If one were to 
incorporate this additional freedom on the value of the 
zero-level (within its error bars), the relative x 2 difference 
would decrease. Therefore we do not feel that the above 
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A\ 2 rules out Oi as a tracer of Hi. In the case of Nl, 
however, our results indicate that it cannot be considered 
as a perfect tracer of H i, as the x 2 difference is quite large, 
Ax 2 = 41.8. Varying the zero- level for the N I lines cannot 
account for this large A% 2 , and we cannot find any other 
explanation than an actual difference in behavior between 
Nl and the other neutrals Di, Oi, and Hi. 

Finally, we wish to point out another important consis- 
tency check: whether we require D I or O I to be a tracer of 
Hi, we obtain the same average D/H ratio. This supports 
the robustness of our measurement of the D/H ratio. 

H I D I 




1214.5 1215.0 1215.5 1216.0 1216.5 1217.0 

WAVELENGTH 

Fig. 9. The final fit for the H I and D I lines when all lines 
are taken into account and all parameters free. 



D I 




1215.0 1215.2 1215.4 1215.6 



WAVELENGTH 



1215.8 1216.0 



Fig. 10. The final fit enlarged over the Di line when all 
lines are taken into account and all parameters free. 

We derived the error bar on the D/H ratio by study- 
ing the relative variation of the \ 2 as a function of the 
(unique) D/H ratio, as in section 3.2. This variation is 



shown in Fig. O, where we also plot the 2a error bars 
estimates corresponding to a A\ 2 = 10 in the upper plot, 
together with the variation of the total H I and D I column 
densities in the lower pannel. We finally obtain the follow- 
ing column densities and D/H ratio with an "effective" la 
D/H error bar: 

N(Hi) To tai = 2.4 ±0.1 x 10 18 cm -2 
N(Di) To tai = 2.68 ± 0.05 x 10 13 c m - 2 
(D/H) Tota i = 1.12 ± 0.08 xlO- 5 , 
(the error bars for N(Hi) and N(Di) are derived by 
trial and error; they are thus more tentative estimates but 
consistent with the D/H error bar). 




1.0 1.2 1.4 1. 
(D/H) Tot (xlO 5 ) 

Fig. 11. Upper pannel: relative x 2 difference with respect 
to the best fit plotted vs. the D/H ratio. A unique D/H 
ratio was assumed on the line of sight. A spline was fitted 
to the numerical results, represented by stars. The dashed 
lines give the minimum and maximum value of the D /H 
ratio that correspond to A^ 2 = 10, or, according to our 
conventions, to an effective 2a confidence interval. Lower 
pannel: variation of the total Hi and Di column densities 
vs. the D/H ratio, as indicated. 



Let us now try two further hypotheses on this global 
line fitting. As we discussed in the previous section, com- 
ponent 3 is identified as the LIC, i.e. the interstellar cloud 
in which the Solar System is embedded. This component 
is also the unique component detected by L93, L95, to- 
ward Capella. Clearly, our average D/H ratio disagrees 
with their estimate of the D/H ratio toward Capella, in 
the LIC: (D/H)= 1.6 ± 0.1 x 10" 5 . We now wish to see 
whether our D/H ratio in component 3, the LIC, can be 
reconciled with this latter value, and thus proceed to pro- 
file fitting, using (D/H)= 1.6 x 10~ 5 in the LIC compo- 
nent. As shown in Table [| this solution is indeed com- 
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Table 6. Physical characteristics of the ISM components, numbered from 1 to 3 in the order of increasing radial 
velocity; component 3 is identified as the LIC, and noted as such. Here the accuracy on the velocities are relative. 
Absolute accuracy on the heliocentric velocities are of the order of 0.8km.s~ . Error bars are tentative la estimates, 
obtained by trial and error. 



Component 


1 


2 


LIC 


V (km.s" 1 ) 


8.20 ± 0.20 


13.20 ±0.20 


20.35 ±0.05 


T(K) 


11000. ± 2000. 


3000. ± 2000. 


4000 + 2000 ' 


a (km.s ) 


0.5 ±0.5 


3.0 ±0.5 


2 0+ 5 - 




0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
(D/H) uc (xlO 5 ) 

Fig. 12. Same as Fig. [OJ except the D/H ratio varied is 
the LIC D/H. In the lower pannel are plotted, as indicated, 
the evolution of the LIC D/H, the D/H averaged over 
components 1 and 2, and the total D/H averaged over 
components 1, 2, and LIC. 



patible with our data set, as the relative % 2 difference is 
reasonable, A% 2 = 5.8. Obviously, such a high D/H ra- 
tio is permitted in the LIC because all three components 
on the line of sight to G191-B2B are thoroughly blended 
in the Di line. It also means that estimates of the D/H 
ratio in a particular component in our dataset is mean- 
ingless, as the D/H ratios between the three components 
are strongly correlated. This effect is shown in Figure |l2| , 
where we plotted the A\ 2 , the D/H ratio averaged over 
components 1 and 2, and the total D/H vs. the LIC D/H. 
Indeed, the error bars on the LIC D/H ratio, that corre- 
spond (at 2cr, according to our conventions) to A\ 2 = 10, 
are rather large. However, it is quite interesting to note 
that the D/H ratio of the LIC, and that averaged over 
components 1 and 2, are anti-correlated. This translates 
itself in that the total D/H ratio remains constant when 
the LIC D/H ratio is varied from 0.6 x 10~ 5 to 2.0 x 10" 5 . 
We obtain the following column densities in the LIC along 



with the corresponding D/H ratio with an "effective" ler 

D/H error bar: 

N(H i) L ic = 1.0 ± 0.3 x 10 18 cm~ 2 
N(Di)lic = 1.1 ±0.1 x 10 13 cur 2 



(D/H) 



LIC 



1.19 



-0.35 
-0.25 



X 10" 



(again, the error bars for N(H i) and N(D i) are derived 
by trial and error; they are thus more tentative estimates 
but consistent with the D/H error bar). 

Finally, if we try to impose all characteristics of the 
LIC found by L93, L95, toward Capella, in our solution, 
we obtain A% 2 =25.6, which seemingly rules out this hy- 
pothesis. As we just showed that the LIC D/H ratio was 
compatible between the two line of sights, we are forced 
to conclude that there must be a significant difference in 
the line broadening of the LIC between both lines of sight 
(see however Section 5). 

Let us therefore conclude on this study, as far as Hi 
and D I are concerned. We find that performing a global 
fit on the line of sight gives good agreement with the indi- 
vidual fits performed on each spectral region. The average 
D/H ratio is tightly constrained to D/H=1.12 ± 0.08 x 
10 -5 . We also found that Di and Oi are very reliable 
tracers of Hi within the quality of our data; this is not 
the case for N I. Moreover, whether D I or O I are set as 
tracers of Hi, i.e. a single column density ratio to Hi is 
required on all three components, then the same average 
D/H ratio is obtained. This constitutes a strong consis- 
tency check of our results. Our average D/H ratio is in 
strong disagreement with that obtained by L93, L95 to- 
ward Capella. The D/H ratio in the LIC, which is the 
interstellar cloud common to both lines of sight, can be 
reconciled with the D/H found toward Capella. However, 
this comes at the expense of a significantly smaller D/H 
ratio in our components 1 and 2, of order 0.9 x 10 -5 , while 
the D /H ratio averaged all three components remains the 
same, (D/H) Total = 1.12 x 10~ 5 . 

4-3. Results for each spectral line 

We now briefly comment on the final solutions for the lines 
of Oi, Ni, Sin, and Sim, as obtained from the previous 
global fit of all lines including H I and D i. These results 
are summarized in Table ^. 

Components 1 and 2 are blended in all observed lines, 
and therefore, individual results for cither one of these two 
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Table 7. Column densities derived in the best-fit solution, obtained simultaneously on all lines and all species. The 
third and fifth column, respectively noted as "1+2", and "Total", give quantities averaged over components 1 and 
2, and averaged over all three components, respectively. These are the most reliable quantities, as components 1 and 
2 are blended in all observed lines. Error bars are tentative estimates, obtained by trial and error. Abundances are 
normalized to 10 6 hydrogen atoms. 



Component 


1 


2 




1+2 




LIC 


Total 




N(Ol) (cm" 2 ) 


2.6 x 10 14 


1.5 x 


10 i4 


4.0 x 


10 i4 


2.9 + 0.3 x 10 14 


6.9 + 0.3 


x 10 


N(Oi)/N(Hi) (xlO 6 ) 
N(Ni) (cm" 2 ) 


823 


135 




288 




2871-° 


288 ± 25 




1.1 x 10 13 


4.6 x 


10 12 


1.6 x 


10 13 


6.3 + 0.4 x 10 13 


7.9 + 0.2 


x 10 


N(Ni)/N(Hi) (xlO 6 ) 
N(Siii) (cm" 2 ) 


36 


4 




11 




RO+33 
°°-18 


33 + 2 




1.9 x 10 13 


8.4 x 


10 12 


2.7 x 


10 13 


1.0 + 0.1 x 10 13 


3.7 + 0.1 


x 10 


N(Siii)/N(Hi) (xlO 6 ) 
N(Sim) (cm" 2 ) 


62 


8 


10 12 


20 


10 12 


10±f 


16 + 1 




3.9 x 10 11 


2.2 x 


2.6 x 


< 1.0 x 10 11 


2.6 ± 0.2 


x 10 


N(Sim)/N(Hi) (xlO 6 ) 


1.3 


2.1 




1.9 




< 0.1 


1.1 ± 0.1 




N(Sin+Sim) (cm" 2 ) 


1.9 x 10 13 


1.1 x 


10 13 


3.0 x 


10 13 


1.0 + 0.1 x 10 13 


4.0 ± 0.1 


x 10 


N(Sin+Sim)/N(Hi) (xlO 6 ) 


63 


10 




22 




10±f 


17+1 




N(Sim)/N(Siii) 


0.020 


0.267 




0.096 




< 0.001 


0.070 





O I 

1.2 1 1 1 1 i 1 1 1 1 i 1 1 1 1 r~ 




1302.15 1302.20 1302.25 1302.30 1302.35 1302.40 
WAVELENGTH 

Fig. 13. The final fit for the Oi line when all lines are 
taken into account and all parameters free. 

components cannot be relied upon. We therefore discuss 
quantities averaged over these two components, and refer 
to them as '1+2'. The LIC characteristics are obtained 
with reasonable accuracy, while total column densities av- 
eraged over the line of sight are always the most accurate. 

Oi 1302 A: 

The fit is shown in Fig. [l3| The reduced \ 2 is 0-88 for 
144 degrees of freedom, i.e. the fit is excellent. 

The total abundance of O I (normalized to 10 6 hydro- 
gen atoms) is N(Oi)/N(Hi)=288±25, almost identical to 
that found for the LIC. This value is also compatible with 
the average ISM gas neutral oxygen abundance, given by 
Meyer et al. (1998a): N(Oi)/N(Hi)=319±14. Meyer et al. 
(1998a) also evaluated that the average ISM (gas+dust) 
should be of the order of 500 assuming all O I atoms have 
returned to the gas phase. For comparison the solar value 



N I 

1.2 1 1 1 1 1 i 1 1 1 1 i 1 1 1 1 r - 




0.0 I i i , i i i , i , i , , i , i , i i , i , i , , i , i , i 

1199.50 1199.55 1199.60 1199.65 1199.70 1199.75 
WAVELENGTH 

Fig. 14. The final fit for one of the N I lines of the triplet 
when all lines are taken into account and all parameters 
free. 

is N(Oi)/N(Hi)=741 ± 130 (Grevesse and Noels, 1993), 
different enough to lead these authors to conclude that 
the solar value may be enriched in oxygen. From Table 0, 
we see that the LIC appears as a normal diffuse interstel- 
lar cloud, as far as O I is concerned, just as the average of 
components 1 and 2. Our observations thus support the 
analysis of Meyer et al. (1998a), with some possible indi- 
cation that component 1 has even returned all its O I to 
the gas phase. 

Nl 1200A, triplet: 

The fit for the strongest line of the triplet is shown in 
Fig. [l4], and the total fit over the triplet gives a reduced 
X 2 per degree of freedom (185) of 1.29. It is certainly less 
satisfactory than the Oi fit; this seems to result from a 
slight mismatch between the three Nl lines. Possibly, this 
may be due to some blending with weak photospheric lines 
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of other species (e.g., Fev, Niv, ...). However, if a fit of 
the triplet is performed independently of all other lines, 
the solution is found to be very similar. We thus remain 
confident in the robustness of the N I results. 

Ferlet (1981) found a strong relationship between 
N(Hi) and N(Ni) in a diffuse Hi medium, on longer 
pathlengths, which translates for 10 6 hydrogen atoms to 
N(Ni)/N(Hi)=62tg. More recently, Meyer et al. (1998b) 
reached the following value for the Ni interstellar abun- 
dance: N(Nl)/N(Hl)=75 ± 4. These evaluations are com- 
patible with the solar value of Grevesse & Noels (1993): 
N(N i) /N(H i)=93 ± 16. We obtain, as an average over the 
line of sight, N(N i) /N(H i)=33 ± 2, in agreement with our 
previous Cycle 1 observations (Lemoine et al. 1996) which 
gave a value of the order of 32, but significatively differ- 
ent from both the solar and ISM values. Interestingly, this 
difference seems to result from components 1 and 2, since 
their average value is ~ 11. Furthermore the LIC nitrogen 
abundance seems to be compatible with the standard ISM 
values: N(Nl)/N(Hl)=63±?|. As we discuss below, com- 
ponents 1 and 2 appear more ionized than the LIC; this 
would then suggest that the value of the Nl/Hi ratio is 
ionization dependent, and could further explain why Ni 
does not turn out to be as reliable a tracer of H I as O I 
and D I were found to be. 

Sin 1304 A: 

The fit is shown in Fig. [l5|. The \ 2 P er degree of free- 
dom (72) is here equal to 0.98 and is thus again highly 
satisfactory. 

We find similar column density ratios in Sin, Mgn, 
and Fen. These two latter spectral regions were observed 
in Cycle 1 at medium and high resolution (Lemoine et 
al. 1996). The correspondence with these previous obser- 
vations is as follows: components 1 and 2 are blended at 
medium and high resolution, and correspond to compo- 
nent A in Lemoine et al. (1996), while, similarly, com- 
ponent 3 corresponds to component B. The detection of a 
component C was reported in these previous observations, 
but it is not confirmed here. It probably was a ghost due 
to the wings of the pre-COSTAR GHRS Echelle-B point 
spread function. As Sin, Fen, and Mgn have similar ion- 
ization properties up to first ionized level, the similarity 
of their column density ratios in each component certainly 
supports our results. We discuss below the ionization and 
depletion of the three components. 

Sim 1206 A: 

The fit is shown in Fig. |l^. The x 2 per degree of free- 
dom (76) is 1.21 and is thus less satisfactory. However, this 
line falls at the end of the spectral range where the S/N 
ratio in the continuum varies strongly from pixel to pixel. 
We did not incorporate this instrumental effect in our pro- 
file fitting procedure. In effect, we adopted a unique S/N 
ratio for the whole continuum (calculated as the average 
over all pixels in the continuum), and derived error bars 



for all pixels using this value of the S/N ratio in the con- 
tinuum, weighted by the flux of the pixel, and combined 
quadratically with the background noise. Therefore, the 
variation of the S/N ratio in the continuum, from pixel 
to pixel, renders our x 2 analysis less meaningful in this 
case. In any case, this line contributes to at most 10% of 
the total degrees of freedom (over all lines), so this cannot 
affect our global fit. 



Si II 




1304.35 1304.40 1304.45 1304.50 1304.55 
WAVELENGTH 

Fig. 15. The final fit for the Sin line when all lines are 
taken into account and all parameters free. 



Si m 




1206.45 1206.50 1206.55 1206.60 
WAVELENGTH 



Fig. 16. The final fit for the Si in line when all lines are 
taken into account and all parameters free. The broad line 
on the red side is the photospheric Si in line observed at 
a 29km. s -1 heliocentric velocity. The LIC contribution to 
this line is negligible. 

We confirm the average N(Sim)/N(Siii)~ 10% found 
in Lemoine et al. (1996). The three components have very 
different N(Sim)/N(Siii) ratios. Components 1 and 2 are 
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more typical of Hn regions, while the LIC is typical of 
an Hi region. The ionization structure of components 1 
and 2 is, however, non trivial. A N(Sim)/N(Siii)~ 0.27 
ratio would indicate a significantly ionized region, yet 
the total broadening (thermal plus turbulent) of the line 
does not confirm this idea. Indeed, component 1 shows a 
high temperature ~ 11000 K and a low microturbulence 
~ 0.5 km.s -1 , whereas component 2 shows the opposite, 
namely, a low temperature ~ 3000 K, and a high microtur- 
bulence ~ 3.0 km.s -1 . These two components may how- 
ever constitute differents parts of a shocked structure, with 
a velocity separation of ~ 5.0 km.s -1 as projected on the 
G191-B2B line of sight. This will be discussed somewhere 
else. 

It is interesting to note that the N(Siin)/N(Sin) ratio 
of the LIC has also been evaluated by Gry & Dupin (1998) 
in the direction of e CMa. They found N(Sim)/N(Sill)~ 

0. 41, a very different result indeed. If confirmed, this differ- 
ence may either show that the ionization structure within 
the LIC is inhomogeneous, or that the LIC may be con- 
fused by coincidence with another ISM component on the 
long e CMa line of sight. More observations are needed to 
clarify this point. 

On the basis of a simple model that calculates the 
ionization equilibrium between Si n and Si in (Dupin and 
Gry 1998), and assuming n(Hll)=n e , we converted the 
Si in / Si II ratios listed in Table 7 into hydrogen ioniza- 
tion ratios, N(Hii)/N(Hi), using the most recent value of 
n(Hl)~ 0.24 cm -3 (Puyoo k BenJaffel 1998). The derived 
value are respectively ~ 0.4, ~ 0.5, and 0.1, for compo- 
nents 1, 2, and the LIC corresponding to n e ~ 0.12, ~ 0.14, 
and 0.03 cm -3 . Apparently, this estimate of the electron 
density of the LIC is much lower than a measurement ob- 
tained toward Capella using the ground and excited state 
lines of Cn (Wood & Linsky 1997): n e = O.llj^cm -3 . 
This discrepancy could be explained by arguing that the 
ionization in the LIC is far out of equilibrium, due for 
example to a shock ionizing event, which would tend to 
enhance the current electron density (Lyu & Bruhweiler 
1996). 

We can also use silicon to probe the degree of deple- 
tion in the three components, as it is more sensitive than 
oxygen or nitrogen. We find a total nitrogen to neutral hy- 
drogen ratio averaged over the three components of ~ 17, 

1. e. a factor two smaller than the cosmic value, ~ 35 (per 
10 6 hydrogen atoms). Therefore, provided the ionization 
of hydrogen is small, as our calculations seem to indicate, 
depletion remains relatively weak. This conclusion can be 
seemingly applied to the three components individually, 
or at least for 1+2 and the LIC, as inspection of Table ^| 
reveals. This low depletion is similar to the one observed 
in hot ISM clouds or in the galactic halo (see e.g. Spitzer 
and Fitzpatrick, 1993, 1995; Fitzpatrick and Spitzer, 1994, 
1997; Jenkins and Wallerstein, 1996; Savage and Sembach, 
1996) as well as in the so called "CMa tunnel" (Gry et ai, 
1995; Dupin and Gry, 1998). 



The comparison between the measured radial veloci- 
ties of the components and the systemic radial velocity 
of G191-B2B, «Gi9i-B2B = 5 ± 2km. s , shows that none 
of these clouds seems to be associated with the star. In 
particular, the association of any of these clouds with an 
expansion shell centered on the star is precluded as the 
clouds have larger radial velocities than G191-B2B. A de- 
tailed discussion of the physical properties is clearly be- 
yond the scope of the present work. This will be the sub- 
ject of a forthcoming paper. 

We derive a total Hi column density of 
N(Hi)=2.4±0.1 x 10 18 cm -2 , and confirm the pre- 
vious evaluation of Lemoine et al. (1996). Our value is 
in disagreement with the EUV measurements of Kimble 
et al. (1993) using the HUT, N(H l)=1.7x 10 18 cm -2 
, Green et al. (1990) using a rocket-borne spectro- 
graph, N(Hi)=1.6xl0 18 cm -2 , Holberg et al. (1990), 
N(Hi)=1.0xl0 18 cm -2 , Bruhweiler & Kondo (1982) 
using Ni and the relationship between N(Nl) and N(Hi), 
N(Hi)=0.8xl0 18 cm -2 . Our estimate is however in 
agreement within the error bars with those of Jelinsky 
et al. (1988) and Paerels & Heise (1989). Moreover, it 
is in good agreement with the recent EUV measure- 
ments of Dupuis et al. (1995) and Lanz et al. (1996), 
N(Hi)=2.1 x 10 18 cm -2 , that use more modern white 
dwarf atmosphere models. 

The fit of the Lyman limit by Kimble et al. (1993), Hol- 
berg et al. (1990) comes through the modeling of the atmo- 
sphere of G191-B2B with homogeneous atmosphere mod- 
els. Their analysis might thus be questioned in view of the 
recent detection of non-negligible amounts of highly ion- 
ized metallic species (Nv, Fev, Siiv) in this photosphere 
(see Vidal-Madjar et al. 1994), although trace of metal- 
lic species are not believed to contribute significantly to 
the absorption over 500A. We find a total Ni column den- 
sity in perfect agreement with that found by Bruhweiler 
& Kondo (1982), and the difference between the N(Hi) 
values is directly related to the difference found in the 
relationship N(Ni)-N(Hi) between our observations and 
that of Ferlet (1981) on different lines of sight. 

5. Discussion 

Our results lead us to the following statement: either 
the D/H ratio is constant in the local ISM, its value is 
1.12±0.08 x 10 -5 , and the estimate made in the direction 
of Capella by L95 is incorrect; or, the D/H ratio does vary 
in the local ISM. If this latter hypothesis is verified, we 
need to check whether the D/H ratio can vary within a 
same cloud, the LIC, which is common to both the Capella 
and the G191-B2B line of sight, or whether it varies on 
larger spatial scales, from cloud to cloud. 

It thus seems reasonable to re-analyze the evaluation 
made by L93, L95 with the same techniques that we used 
to derive the D/H ratio toward G191-B2B. This we do 
now. 
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5.1. Capella revisited 

The Capella high resolution GHRS Echelle-A data are of 
high quality, with a spectral resolution AA = 3.5km.s _1 , 
and a signal-to- noise ratio S/N~ 40 at the interstellar Ly- 
man a line continuum. For this star, the accuracy on the 
estimate of the D/H ratio is limited by the ability to re- 
construct the Lyman a stellar emission profile of Capella, 
all the more since it is composed of the combined profiles 
of the two cool stars of the binary system. L95 circum- 
vented this problem by observing Capella at two different 
phases, and were thus able to reconstruct, at least partly, 
both stellar continua. 

As our paper is not dedicated to the study of this line 
of sight, we only summarize our study of the Capella data 
set, and jump to the conclusions. In order to model the 
stellar continuum, we used two different approaches, in the 
same spirit as Section 3.1: in the first approach, we used 
the stellar continuum produced by L95, and in the other, 
we constructed a stellar profile in a very naive manner, by 
interpolating the far wings of the emission profile with a 
7th degree polynomial. This latter polynomial would then 
be further corrected in the core by a 4 th order polyno- 
mial. In all our fits, we found that the continuum found 
by L95 was the most appropriate in terms of \ 2 statis- 
tics, and hence we will not refer to this 7 th x 4 th order 
polynomial continuum anymore. L95 detected only one 
component on the line of sight, from the analysis of the 
narrow metal lines of Mgn and Fen. We followed these 
authors and proceeded to the profile fitting with only 
one component on the line of sight. A first analysis gave 
results in agreement with their work, although we find 
N(Di)= 2.6 ± 0.1 x 10 13 cm~ 2 , in disagreement with the 
original estimate of L95, N(D i)= 2.85 ± 0.1 x 10 13 cm~ 2 . 
It probably results from a slight mismatch in the radial 
velocity of the Di and Hi lines, with a shift of 0.3km.s _1 
for the D I line away from the expected position, to be 
compared with an error on the radial velocity of only 
i0.lkm.s~ 1 . Whereas our Di and Hi lines automatically 
have the same radial velocity, L95 did not link the radial 
velocities of the D I and of the H I absorbers. Nevertheless, 
this seems to point to the existence of another component, 
that would perturb the H I line but not the D I line. 

Wood & Linsky (1998) did mention that such a per- 
turbation should be present, and due to the absorption by 
the hydrogen atoms in the Earth geocorona, which, at the 
time of the observations, was located in velocity space a 
— 22km. s _1 . This absorption is shown in Fig. [l7]. When we 
attempt a fit with a two component solution, we obtain a 
better x 2 with a A\ 2 ~ 100, a dramatic improvement in- 
deed. In this two component solution, we left both compo- 
nents entirely free. Interestingly enough, we obtained two 
kinds of solution: in both cases, the main component was 
extremely similar to the single component (LIC) of L95; 
in one case, however, the extra absorption was caused by 
a cold, i.e. small broadening component, located around 



Table 8. Physical parameters evaluated on the line of 
sight to Capella from Cycle 1 GHRS archival data. We 
assumed the presence of three components on the line of 
sight, which represents a gain A\ 2 ~ 100 in \ 2 over the 
one cloud solution. The additional components are noted 
"Hot" and "Geo" and were left as entirely free in the pro- 
file fitting. The data are here normalized with the stel- 
lar continuum reconstructed by L95 and corrected by a 
4 th order polynomial. 



Stellar continuum 


Linsky x4 ttL order polynomial 


7? /198 


0.79 


N(Hi)lic (cm" 2 ) 


1.72 x 10 18 


N(Dl) L ic (cm" 2 ) 


2.69 x 10 13 


(D/H) L ic (xlO 5 ) 


1.56 


Vlic (km.s- 1 ) 


23.05 


T L ic (K) 


6053. 


ctlic (km.s -1 ) 


2.96 


N(Hl) Gco (cm" 2 ) 


1.74 x 10 12 


Vgco (kirns -1 ) 


-22.00 


Tgco (K) 


2113. 


N(Hl) Hot (cm" 2 ) 


4.00 x 10 14 


V Ho t (km.s -1 ) 


+21.50 


T Ho t (K) 


26667. 



~ — 20km. s ; in a second case, slightly favored over the 
previous one, in terms of x 2 , we obtained a hot component, 
with T ~ 30000K, located around ~ +20km.s -1 . Our pro- 
file fitting algorithm indeed works by random minimiza- 
tion, and can therefore find different degenerate solutions 
for a same dataset. Since both of these extra components 
fall at smaller radial velocities than the LIC, they produce 
a shift of the H I absorption, and reconcile the overall ra- 
dial velocities of both H I and D I lines. 

Since we know that the geocorona has to be present, 
we performed a three cloud solution, with all parameters 
free. We found, in that case, both previous extra compo- 
nents, i.e. one hot, located around ~ +20km.s~ , and one 
cold component, located at — 22. 0km. s -1 ! The final char- 
acteristics are given in Table g and the fit is shown in 
Fig. 0. The geocorona found here is entirely consistent 
with the actual geocorona, in terms of H I column density, 
velocity, and temperature. 

The "Hot" component falls at an heliocentric veloc- 
ity of order 21.5 ± 1.5km.s — , with an effective broaden- 
ing temperature T ~ 27000i4QQQ°K. Its characteristics are 
well defined, even though it seems to be constrained from 
one side only of the H I line, probably because its pres- 
ence is felt, in these high quality data, over nearly lA. 
This component may be identified with a cloud interface, 
as first detected by Bertin et al. (1995) toward Sirius A, 
or with a "hydrogen wall" as proposed and discussed by 
Linsky & Wood (1996), Piskunov et al. (1997), Dring et 
al. (1997), and Wood & Linsky (1998). 

Its slightly shifted value redward of the LIC is com- 
patible with the "hydrogen wall" concept since it assumes 
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CAPELLA (ISM + Geocorona + Hot) 
1.2 I 




1215.40 1215.50 1215.60 1215.70 1215.80 1215.90 1216.00 
Wavelength 

Fig. 17. The best fit obtained with three components as- 
sumed along the Capella line of sight. The stellar contin- 
uum assumed here is the one reconstructed by L95 and 
corrected by a 4 th order polynomial. The two additional 
components are imposed by the data points near 1215.6A 
which cannot be fitted by any simpler "one cloud" solu- 
tion. The weak additional component needed at -22km. s -1 
is at the velocity shift of the earth geocorona but cannot 
alone fit correctly the 1215. 6A feature. A third "hot" com- 
ponent is also needed (see text). 



a slowing down of the ISM hydrogen atoms when inter- 
acting with the stellar wind itself. If correct, this seems to 
indicate that the LIC extends out to Capella, and that it 
is quite homogeneous in that direction, since the average 
volume density toward Capella is then very similar to the 
one observed within the solar system. 

We expect this component to be entirely compatible 
with our observations of G19TB2B as it would fall in 
between components 2 and 3, in a region where the Hi 
Lyman a line is completely saturated. 

In summary, our estimate of the D/H ratio toward 
Capella is: 

(D/H) CapcU a = 1-56 ± 0.1 xlO- 5 , 
in perfect agreement with the original estimate. From 
these data alone we also estimated the other characteris- 
tics of the LIC independently and found that its tempera- 
ture should be Tlic = 6750±250 K if the microturbulence 
of the LIC is properly constrained by the other spectral 
lines as given in L95. 

As a conclusion, we thus confirm the estimate of L95 of 
the D/H ratio within the LIC, (D/H)= 1.6±0.1 x 1CT 5 , but 
seem to evaluate a slightly cooler LIC temperature, more 
compatible with the G191-B2B evaluation. The LIC D/H 
estimate made toward Capella is thus marginally com- 
patible with the evaluation made within the LIC in the 
direction of G191-B2B, and in any case clearly incompat- 



ible with the estimate of the average D/H ratio on this 
latter line of sight. 

The D/H ratio thus has to vary within the local ISM. 
6. Conclusion 

The outcome of our Cycle 5 observations of the white 
dwarf G191-B2B performed with the GHRS of HST, aim- 
ing to derive "bona fide" value(s) of the D abundance in 
the local interstellar medium is the following : Overall 8 
absorption complexes of the elements Sin, Sim, Ni, Oi, 
and, of course, Hi and Di, were observed at high reso- 
lution (Echclle-A grating), and analyzed simultaneously 
using a new profile fitting procedure. We extract an ac- 
curate description of the cloud structure along the line of 
sight, which consists of three absorbing clouds, separated 
by ^5 km.s -1 . The third absorber, to be understood in 
terms of increasing radial velocity, at a heliocentric ve- 
locity V — 20 .4 ± 0.8km.s _1 , is identified with the Local 
Interstellar Cloud in which the Sun is embedded. The net- 
work of constraints on these clouds, that results from the 
profile fitting of the different interstellar lines, is such that 
we can derive the temperature and the turbulent velocities 
in the different interstellar components, as well as the col- 
umn densities of the various observed species. Our new set 
of high quality data is consistent with the Lemoine et al. 
(1996) evaluation of the total Hi column density towards 
G191-B2B, i.e. N(Hl)=2.4±0.1 x 10 18 c m - 2 . 

We derive an average D/H ratio over the three ab- 
sorbing clouds N(Di) total /N(Hi) tota i=1.12±0.08 x 10- 5 . 
The Local Interstellar Cloud, detected here toward G191- 
B2B, has also been detected on the line of sight toward 
Capella, and its D/H ratio has been measured by Lin- 
sky et al. (1993, 1995): (D/H) CapcUa = 1.6 ± 0.1 x 10" 5 . 
We have re-analyzed the data of L95 toward Capella, 
and confirm their estimate of the D/H ratio, as we find: 
(D/H) Capclla = 1.56 ± 0.1 x 10~ 5 . We find that the D/H 
ratio in this Local Interstellar Cloud, on the line of sight to 
G191-B2B can be made consistent with the Capella value. 
However, this comes at the expense of a much smaller av- 
erage D/H ratio in components 1 and 2, of order 0.9 x 10 -5 , 
in such a way that the D/H ratio averaged over all three 
components remains at the above value 1.12 x 10~ 5 . 

Therefore, we conclude that the D/H ratio indeed 
varies in the local ISM over a few parsecs and at least 
from interstellar cloud to interstellar cloud. Although we 
do not detect variations of the (D/H) ratio within the Lo- 
cal Interstellar Cloud itself, such variations may become 
clear as more data are analyzed. Here, we note that both 
clouds 1 and 2 are more like Hn regions while the Local 
Interstellar Cloud could be identified with an Hi region; 
this ionization difference might point to the cause of this 
(D/H) variation. 

These large variations of the ISM D/H ratio in so small 
scales cannot be due to any nuclear processes (spallation 
processes induced by cosmic rays or by huge gamma ray 
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fluxes) occuring in that medium. We mention here two 
other possibilities: 

(i) - given the fact that the column density is only of 
the order of a few 10 18 cm -2 the lowest D/H value could 
be due to the admixture of a comparable amount of D- 
free, pure H remnant of a material within the considered 
cloud. 

(ii) - these low D components could be the remnants of 
semi-thermalized blobs of material coming from the outer 
layers released during the planetary nebula phase of evo- 
lution. 

As a rule, it is much easier to deplete than to create 
deuterium. In this case, just as in the two possibilities 
mentioned above, the lowest ISM D/H values would be 
due to a secondary effect, and the highest value would be 
more representative of an actual ISM D/H ratio. However, 
nothing precludes the Capella value from being itself af- 
fected by such a secondary effect. Moreover, we point out 
that the deuterium abundance may be locally enhanced 
(or depleted) when the radiation flux is anisotropic, as the 
radiation pressure acts differently on H I and D I (Vidal- 
Madjar et al. 1978; Bruston et al. 1981). 

To sum up the consequences of the present investi- 
gation, the question concerning the existence or not of 
variations of the D/H ratio in the ISM has now an an- 
swer : variations do exist within the local ISM of order 
« 30% — 50% within a few parsecs, and/or within a given 
cloud, the Local Interstellar Cloud. To be more precise, 
our observations show that the average (D/H) ratio on 
the line of sight to G191-B2B (averaged over three in- 
terstellar components), is significantly different from the 
(D/H) ratio measured toward Capella. Although the Lo- 
cal Interstellar Cloud is common to both lines of sight, we 
do not find conclusive evidence for variation of the (D/H) 
ratio in this particular component. 

The comparisons between the theoretical predictions of 
standard Big-Bang nucleosynthesis combined with those 
coming from the recent models of chemical evolution of 
deuterium and the measured D/H ratios in the different 
components of these lines of sight appear to be more com- 
plex to establish. Forthcoming analyses will have now to 
take into account these D /H variations and should discuss 
more thoroughly some of the possibilities envisaged above. 

We therefore look forward to: (i) further observations 
of the interstellar deuterium abundance on other nearby 
lines of sight. In this respect, white dwarfs seem to con- 
stitute a highly reliable way of achieving a precise mea- 
surement of the deuterium abundance in the nearby in- 
terstellar medium. In effect, white dwarfs present nei- 
ther the complex stellar continuum of cool stars, nor the 
complex interstellar velocity structure that arises towards 
more distant hot stars, (ii) Measurements of the D/H ra- 
tio in spiral arms of the Galaxy distinct from ours, or in 
low-redshift systems, performed in the framework of the 
FUSE-LYMAN mission. 
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